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Recently it was shown that kink baryons and kink-antikink scattering in the massless Gross-Neveu 
model are closely related to one- and two-soliton solutions of the sinh-Gordon equation. Here we 
generalize these findings to the case of n kinks and antikinks. Using the known n soliton solution of 
the sinh-Gordon equation, we solve the general n kink-antikink scattering problem in the large N 
Gross-Neveu model analytically, mapping the time-dependent Hartree-Fock approach onto inverse 
scattering theory. 

PACS numbers: 11.10.-z,11.10.Kk 



I. INTRODUCTION 

In this paper, we continue our study of the simplest 
Gross-Neveu (GN) model [![, a 1+1 dimensional model 
field theory of N species of massless, self-interacting 
Dirac fermions with Lagrangian 

N 2 / N \ 2 

fc=i / 

We restrict ourselves from the outset to the 't Hooft limit 
N — > oo, Ng 2 = const. Semiclassical methods have re- 
vealed a number of fascinating properties of this model 
over the years, see the review articles [2-4] and references 
therein. A key quantity in these studies is the scalar 
mean field S. It plays a role similar to Witten's "mas- 
ter field" in large N gauge theories [f|, namely as saddle 
point of the functional integral from which all observables 
can be computed. For fermions in the large N limit, it 
can be identified with the self-consistent Hartree-Fock 
(HF) potential. 

Most of the results for S obtained so far are related to 
static problems. In the vacuum, the HF potential is ho- 
mogeneous and can be interpreted as dynamical fermion 
mass [l|. Localized, spatially varying HF potentials de- 
scribe individual baryons [(| . Spatially periodic solutions 
appear in investigations of baryonic matter, both at zero 
and finite temperature Q . The most difficult problem 
is to find solutions of the time-dependent Hartree-Fock 
approach (TDHF), at least non-trivial solutions which 
are not simply boosted, static solutions. The only known 
analytical solutions of this type to date are the breather 
Q and kink-antikink scattering Q. Since both arc re- 
lated by analytical continuation, there is in fact only one 
non-trivial time-dependent solution known. This reflects 
the lack of systematic methods to derive time-dependent, 
self-consistent mean fields for fermions. 

Recently, it was pointed out that the situation is more 
favorable for a class of particularly simple TDHF solu- 



* fitzncr@thcorie3 . physik. uni-crlangon. dc 
' thics@thcoric3 . physik. uni-crlangcn.dc 



tions, classified as "type I" in They are defined as 
those solutions where the scalar density of each single 
particle level is proportional to the full self-consistent po- 
tential S, 

4>a4>a = A a 5, (2) 

where A Q may vanish for some states. If property ([2]) 
is satisfied, the TDHF problem reduces to the classi- 
cal N = 1 GN model, for which Neveu and Papanico- 
laou have uncovered a relationship with the sinh-Gordon 
equation some time ago [lol ]. As a consequence, the self- 
consistent TDHF potential of the GN model ([1} can be 
shown to satisfy the classical sinh-Gordon equation 
This is surprising at first sight, as the sinh-Gordon equa- 
tion possesses only singular solitons. Owing to a non- 
linear field transformation however, these singularities 
are mapped onto zeros of S, 

□it + 4 sinhu = 0, u = lnS 2 , (3) 

so that the scalar mean field S is perfectly regular. One 
can easily check that the mean fields for the kink baryon 
PH , kink-antikink scattering 0, Q and the kink crystal, 
the ground state of the GN model at finite density [3], 
are indeed all related to known soliton solutions of the 
sinh-Gordon equation. 

This raises immediately the question: Are there other 
soliton solutions of the sinh-Gordon equation which 
might yield physically sensible, new TDHF solutions of 
the GN model? If one thinks about this problem, one en- 
counters two potential obstacles. The first has to do with 
the singularities of all sinh-Gordon solitons, the second 
with the fact that the sinh-Gordon equation is a nec- 
essary condition for type I solutions, but perhaps not 
sufficient. 

The first difficulty can be handled as follows. If one 
inspects the available solutions of the sinh-Gordon equa- 
tion in the literature, one finds in all cases that the ar- 
gument of the logarithm in Eq. ([3]) has either zeros, or 
poles, or both. This reflects the fact that all solitons 
are singular. In order to get a bounded S, we should 
only allow for zeros. As already pointed out in the 
most interesting solution of this type is presumably the 
n soliton solution constructed by the inverse scattering 
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method [12hT4T ] (throughout this paper we use N to de- 
note the number of flavors and n to denote the number 
of baryons, to avoid confusion). In the GN model the n 
soliton solution is expected to describe time dependent 
scattering of n kink- and antikink-baryons. 

The second difficulty simply means that solutions of 
the sinh-Gordon equation should only be taken as candi- 
dates for TDHF potentials in the GN model. Given any 
such solution, one still has to solve the Dirac equation 
for all continuum states in the Dirac sea and the valence 
bound states and verify self-consistency of the mean field. 

In this paper, we propose to do just this for the n soli- 
ton solution. Our main goal is to find the most general 
type I solution of the TDHF equations for the GN model. 
From the particle physics point of view, one is rarely in- 
terested in scattering problems involving more than two 
incident particles. A time-dependent n soliton solution 
on the other hand describes a scattering process involving 
n incident and n outgoing objects. As a purely theoret- 
ical problem, we find it nevertheless challenging to solve 
the dynamics of n composite, relativistic bound states at 
the elementary fermion level, in full generality. Our mo- 
tivation is not primarily particle physics phenomenology, 
but the desire to find new exact, analytical solutions of 
a relevant model quantum field theory. 

Finally, let us try to relate our work to another impor- 
tant property of the GN model, integrability. As is well 
known, the GN model (TTJ is an example of an integrable 
quantum field theory for any value of N. The exact S ma- 
trix, including kinks and antikinks, has been constructed 
some time ago Nevertheless we find it worthwhile 

to attack this problem with entirely different methods in 
the large N limit. First of all, the S matrix for the finite 
N GN model is only known in principle. The examples 
worked out in the original references deal with low values 
of N (2-8 Majorana flavors, corresponding to 1-4 Dirac 
flavors) and few particles only. Since the algebraic com- 
plexity rapidly increases with increasing number of fla- 
vors and participants, it is not easy to infer the large N 
limit of the collision of n bound states from the published 
S matrix. Secondly, the full dynamical TDHF solution 
has more information than the S matrix which encodes 
only asymptotic, on-shell scattering information. Finally, 
although integrability certainly helps to find the TDHF 
solution, it is apparently not a prerequisite. Thus for in- 
stance, although the massive version of the GN model is 
not integrable, HF solutions have been found for baryons 
[HI, [lj| and baryonic crystals 0, [2l[ in closed analytical 
form. For all of these reasons we have decided to make 
a dedicated effort to solve the n kink-antikink scattering 
problem in the large N limit of the GN model. 

The paper is organized as follows. In Sec. II, we give a 
rather detailed introduction into the single kink baryon 
in an arbitrary Lorentz frame and set up our notation 
in light-cone coordinates. Sec. Ill briefly recalls the n 
soliton solution of the sinh-Gordon equation. In Sec. IV 
we describe how we get to the TDHF spinors and prove 
self-consistency. Sec. V is needed to put the formal re- 



sults into a form better suited for practical applications, 
which then follow in Sec. VI. Here we characterize the 
general n baryon scattering process qualitatively and ex- 
hibit a few illustrative examples involving dynamics of 
up to eight solitons. We end with a concluding section. 

II. REVIEW OF THE KINK BARYON 

The kink baryon of the GN model, originally discov- 
ered by Callen, Coleman, Gross and Zee (cited in [ll| ) , is 
at the same time the simplest and the most exotic baryon. 
Its properties are well studied [f| [22| - [25j . We review it 
here because of its role as main actor in the dynamical 
n baryon problem addressed in this work. An important 
aspect in which we differ from all previous works except 
[25[ is the fact that we consider the kink in an arbitrary 
Lorentz frame, not just its rest frame. This is of course 
a prerequisite for treating scattering processes. 

The outline of this chapter is the following: We will in- 
troduce light-cone coordinates and present first the vac- 
uum and then the boosted kink in the TDHF approach. 
The scalar HF potential S and the self-consistency issue 
will be addressed. We then compute expectation values 
of other relevant fermion bilinears, namely the density 
p = ip'ip, the pseudoscalar density tp'i'j^tp and the ax- 
ial charge density ps = ip^j^ip, resolving contributions 
from the Dirac sea and the bound state. Next we briefly 
recall the derivation of the sinh-Gordon equation from 
Ref. for type I TDHF solutions, of which the kink is 
a paradigm. Finally we summarize the essential physics 
properties of the kink. This section presents no new re- 
sults, but serves to introduce light-cone coordinates and 
set up the notation to be used in later chapters for the n 
baryon problem. 

Starting point is the TDHF equation of the GN model, 
expected to become exact in the large N limit, 

occ 

(z 7 M 9 M -S)ip Q = 0, S = -g 2 M a - (4) 

The sum over occupied states runs over the whole Dirac 
sea as well as possible valence states and includes fla- 
vor degrees of freedom. A non-vanishing scalar mean 
field S signals breakdown of the Z2 chiral symmetry tjj — > 
75^, ipi/) —4>ip- We choose a chiral basis for the Dirac 
matrices, 7 = oi,^ 1 — 102, where 75 = 7°7 1 = — 03 is 
diagonal. In conjunction with light cone coordinates 

z = x — t, z = x + t, do = 8 — d, d\ = 8 + d, (5) 

this simplifies the Dirac equation in (j¥]) to 

2i8ip 2 = Sipi, 2idip! = -Sip 2 (6) 

in terms of upper, left-handed (ipi) and lower, right- 
handed (^2) spinor components. 

Consider first the vacuum problem where S = m = 1 
is the dynamical fermion mass in natural units. Here, 
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the TDHF equation reduces to the free, massive Dirac 
equation with solutions 



A = 



1 



(7) 



labeled by a spectral parameter £. This parameter 
contains the information on momentum k and energy 
u> = ±\/fc 2 + 1 via 



k = C 1 oj = —C , 



(8) 



a relation which allows us to cast the plane wave factor 
in (|7|) into the standard form 



a i(kx— ujt) 



(9) 



(If z is interpreted as light cone time, then £ is the 
light cone energy, but we shall not use this language in 
the following.) The gap equation arises from the self- 
consistency equation for the scalar condensate in the vac- 
uum. The continuum spinor tp^ yields the scalar density 



tp c tp c = -- 



4c 



4C 2 



(10) 



The (cutoff regularized) summation over the Dirac sea 
can be performed conveniently after the following change 
of integration variables, 



A / 2 dk 



A/2 



2 71 



A/2 d( - 1 + 4C 2 
1/2A2TT 4C 2 



The resulting gap equation 

2 f A/2 d 
1 = Ng 2 1 



Ng 2 



In A, 



(11) 



(12) 



1/2A 27r C 

yields the relation between bare coupling and cutoff char- 
acteristic for dimensional transmutation. 

We now turn to the simplest baryon solution of 
Eqs. j4|, the kink or antikink. Without loss of gener- 
ality, we consider the antikink moving with velocity v\. 
In ordinary coordinates it is given by 

S=- tanh (7x(x - v x t) + a x ) , ^ = 1 (13) 

V i - v i 

interpolating between the vacua S = 1 at x — > — oo and 
S = — 1 at x —> +oo (the results for the kink —S can 
simply be generated by a 75 transformation). 

In what follows, it will be advantageous to express 5* 
through exponentials, 



S 



1-Tl 
1+Tl 



ti = exp{27i(x - vit) + 2a{\ . (14) 



Switching to lightcone coordinates, the basic building 
block, n, can be seen to be closely related to a "plane 
wave" with imaginary spectral parameter, 



= p l ( <lZ -3fr) 



+a 1 



-7l(l - Vl). 



(15) 



This structural clement will be important later on. The 
TDHF spinors for the antikink can easily be found. In 
lightcone notation, the continuum states read 



1 



v/rwv + ) i+n ' (b) 



where k\ and t\ differ only by a constant, complex phase, 



Ci-C 
Ci + C 



71. 



(17) 



As is well known, the potential S is reflectionlcss, a cru- 
cial property for everything we shall do in this work. 
The kink at rest possesses one normalizable zero energy 
bound state, in agreement with the expectation based 
on its topological properties. The corresponding boosted 
bound state can be obtained from the continuum spinor 
by setting C — Ci (i- e -i analytic continuation to imaginary 
spectral parameter) and normalizing, 



1 



2& 
-1 



(18) 



/2zCT V -1 / i + n 

The scalar densities for continuum and bound states, 
4C - 



1 + 4C 2 



0. 



(19) 



show that we are dealing with a type I solution accord- 
ing to the classification of Ref @ — every occupied state 
yields a contribution to the scalar condensate propor- 
tional to the full HF potential S. The self-consistency 
condition simply reduces to the vacuum gap equation 

m, 



A/2 



1/2A 



d( 1 + 4C ; 
2tt 4C 2 



^InA. 



(20) 



Consider the expectation value of the fermion density 
in the kink next. It consists of two contributions, one 
from the continuum states (the Dirac sea) and one from 
the bound state. An individual continuum state ip^ gives 
the following (vacuum subtracted) contribution to the 
density 



V^c - 1 



4C 2 (1-4C 2 ) 



(l+4C 2 )(d 2 -C 2 )(l + ri) 2 



(21) 



Performing the d^-integration and multiplying by the 
number of flavors (each state is fully occupied), we find 
the continuum fermion density 



Pcont 



N 



-JV71 



1 dC 1 + 4C 2 
2tt 4C 2 

(l + n) 2 = 



«M - 1) 
Nd x S 



(22) 



and hence the following contribution from the Dirac sea 
to the total fermion number, 



A -JL 
(XX Pcont " 2 ' 



(23) 
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This result can be understood hcuristically as follows: 
The midgap state receives one half of its strength from 
the negative, the other half from the positive energy con- 
tinuum. This half missing state in the Dirac sea manifests 
itself in the peculiar value of the induced fermion num- 
ber ([23)1 . This effect has been discussed extensively in 
the context of fractional fermion number and gives rise 
to observable consequences in condensed matter systems, 
such as unusual spin-charge assignments in solitonic ex- 
citations of polymers [26l - l28{ . 

Next we turn to the contribution to the fermion density 
from the bound state, assuming that the valence level is 
filled with Ni(< N) fermions. The bound state fermion 
density is 



(i + n) 2 ' 



(24) 



normalized to the number of fermions in the valence 
state. 



dxp 



(l) 



(25) 



The continuum and bound state densities (|22[) and (pM)) 
are proportional to each other, so that the total fermion 
density becomes 



(p) = Pcont + = ( 



The total fermion number Ni — N/2 of the kink lies be- 
tween —N/2 and +N/2. In particular, if the bound state 
is half filled, the density vanishes identically. We are then 
dealing with a time-dependent excitation of the scalar 
condensate, a pure "domain wall" moving with constant 
velocity V\ . If the bound state is fully occupied or empty, 
the kink carries ±N/2 fermions and may be thought of 
somewhat loosely as half a baryon or antibaryon. 

For the sake of completeness, let us also evaluate the 
pscudoscalar condensate along similar lines, once again 
assuming N\ valence fermions, 



(^75^) = ^ 



N 
~2 



2n 



(l + n) 2 



(27) 



This quantity is finite and vanishes in the vacuum, so 
that no subtraction is needed. Finally, the last inde- 
pendent bilinear is the axial density (or vector current) 
ps = ip'j5ip, where we must once again subtract the vac- 
uum contribution, 



(P6> = ( N! - ^ ) 2j!Vi 



(28) 



Notice that in all 3 cases (|26l27l28p . the sum over con- 
tinuum states is proportional to the contribution from 
the bound state, with identical relative weights (discrete 
and continuum parts can be identified via the factors iVi 
and N, respectively). This fact can be understood with 



the help of the divergence of vector and axial vector cur- 
rents 12911 . 



Invoking large ./V factorization and using 



3° = 35 = P> 3 1 = 35 = Pz, 



characteristic for 1+1 dimensions, we get 



do(p) 
d (p5) 



9i (pb) 



0, 

25(^75^), 



(29) 



(30) 



(31) 



showing that the three bilinears p, p^^ipij^ip are linearly 
related. As a test of the above calculations, one can verify 
that the kink results for the bilinears do satisfy Eqs. pip . 

The evaluation of mass, energy and momentum of 
the kink baryon is delicate due to vacuum subtraction 
and subtleties in the counting of modes. We refer to 
Ref. (2f| where it was shown in detail that the TDHF ap- 
proach gives a covariant energy-momentum relation for 
the baryon in the GN model, 



M = — , 

7T 



E = 7i M, P = 71U1M 



(32) 



(in natural units). The mass of the kink is independent 
of the number of fermions carried by it, since the bound 
state has zero energy in the rest frame and vanishing 
chiral condensate. 

So far, we have only dealt with the Dirac equation 
involving Btp2 and dipi- As shown in the other two 

derivatives, dipi and dtp2, can also be expressed linearly 
in ipi , tp2 with coefficients depending on S and its first 
derivatives. The result, valid for type I solutions if S 
approaches a vacuum value ±1 for x — > ±oo, is a kind of 
"extended Dirac equation" 



dip = dip, dip = C 2 ip 



(33) 



with 



Ci = 



c 2 = 



S^dS -2iC 2 S^ 
-iS/2 

iS/2 

s^s- 1 s- l as 



The integrability condition of the system 



dC x -dC 2 + [Ci,C 3 ) - 
yields the sinh-Gordon equation for u ■ 

ddu — sinh u = 
or, in normal coordinates, 

□u + 4 sinh u = 0. 



0. 

= \nS 2 



(34) 



(35) 



(36) 



(37) 



The linearized form of this last equation is the Klein- 
Gordon equation for a scalar field with mass 2 which 
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may be identified with the well-known a meson of the 
GN model. Hence the kink can be thought of as a large 
amplitude excitation of the a field, thereby extending the 
Skyrmc picture to the case of a discrete chiral symmetry. 
Finally, it is easy to check that Eqs. (|33I34|) hold for 
continuum states (real £) as well as for the bound state 
(imaginary C, C = &)• 

Summarizing, let us enumerate some properties of the 
kink which will turn out to be important for the case of 
n interacting kinks as well: 

1. The TDHF solution is reflectionless and of type I. 

2. There is a single bound state with vanishing scalar 
density, related to the continuum states by analytic 
continuation in the spectral parameter. 

3. The contributions to the fermion density from the 
continuum states and the bound state have the 
same functional form. The fermion density van- 
ishes identically for a half filled valence level. 

4. Shape, mass and motion of the kink are indepen- 
dent of the number of fermions it carries — in this 
sense, there is no backreaction of the fermions. 

We should like to point out that in spite of the solvabil- 
ity of the model and the peculiar properties of the kink, 
we are dealing with a rclativistic, composite object with 
an interesting internal structure reminiscent of hadrons. 
In Ref. [25|, the structure function, derived analytically 
from the fermion momentum distribution in the infinite 
momentum frame, was shown to display non-trivial con- 
tributions from "valence quarks" , "sea quarks" and "an- 
tiquarks" , with a slight abuse of language. It is therefore 
a non-trivial question to ask how such composite, rela- 
tivistic objects interact with each other. 



III. MULTI-SOLITON SOLUTION OF THE 
SINH-GORDON EQUATION 

As discussed above, the kink of the GN model is akin 
to the one-soliton solution of the sinh-Gordon equation. 
Similarly, the kink-antikink scattering problem can be 
mapped onto the two-soliton solution [9j . If the n baryon 
TDHF solution S of the GN model is of type I, then In S 2 
must also be a solitonic solution of the sinh-Gordon equa- 
tion. An obvious candidate is the known n soliton solu- 
tion of the sinh-Gordon equation, constructed with in- 
verse scattering methods [L2h14| ] . Here we collect all for- 
mulae needed to solve the n baryon problem later on. We 
closely follow the notation of Jevicki and Jin [l4[ . Since 
the focus of our work is not on classical soliton theory 
itself but rather on the role solitons play in the TDHF 
approach, we postpone the discussion of the physics to 
Sec. EII1 

It is inherent in the inverse scattering method that the 
soliton solution of a nonlinear partial differential equa- 
tion is accompanied by a linear problem involving 2- 
componcnt "spinors" . These auxiliary spinors depend on 



a spectral parameter £. In the case of the sinh-Gordon 
equation, they are given in light-cone coordinates (|5|) by 



^(C,z,z) = -A T (0 
<P2(C,z,z) = (l + A T (C)a 



1 



I -a 2 
1 



X e 



I — a 2 

Here, A and A are n component vectors, 

X k = ./TTTnV^ 2 -^) 
A*(C) 



(38) 



7^(0)e 
Xk 



whereas a is the symmetric n x n matrix 

XkXi 



aid 



(39) 



(40) 



Cfc + 

The spinor ip satisfies the system of differential equations 
dip = Uip, dip = Vip (41) 

with 



— i( Tjdu 



V 



i I cosh u — sinh u 
4£ I sinh u — cosh u 



(42) 



u is the solution of the sinh-Gordon equation 

ddu — sinh it = 0, (43) 

as can be shown with the help of the integrability condi- 
tion 



dU - dV + [U, V}=0 
and is related to tp via 

'4Cd(ip! + ip2) 



= In 



(44) 



(45) 



It does not depend on the spectral parameter £, as can 
be seen more easily from the equivalent expression 



In 



det 



I- a 
1 + a 



n 2 



(46) 



Like all soliton solutions of the sinh-Gordon equation, 
the function u of Eqs. (|45|46[) is singular — in fact the n 
soliton solution has n singularities. We identify e u with 
S 2 , the square of the TDHF potential, and will derive 
the TDHF wave functions from ip\,ipi- In this process, 
singularities of u are mapped onto zeros of S which is 
bounded. By comparing e u with S 2 in the one soliton 
case, we can identify the parameters (fc,Cfc(0) as follows 
[see Eq. (T5j)], 



Ik 



1 



Cfc(0) = 2C fe e 



2a, 



(47) 



G 



Vk is the (asymptotic) velocity of the k-th soliton, ajt 
is related to its initial position. Hence the solution is 
general enough to describe the n soliton problem with 
arbitrary initial positions and velocities of the solitons. 
Furthermore, one can verify that tp±, p 2 satisfy 



= -i 



(48) 



for all n. Indeed, by differentiation one finds that the left- 
hand side is independent of z, z, using Eqs. (|41I42[) . The 
integration constant can be taken from the asymptotic 
region. Property (|48[) will be crucial for the proof of self- 
consistency in the following section. 



IV. TDHF SOLUTION FOR N BARYON 
SCATTERING VIA GAUGE TRANSFORMATION 

The sinh- Gordon equation provides us with candidates 
for the simplest class of TDHF solutions (type I) of the 
large N GN model. In each case one still has to verify 
self-consistency of the result. To this end one has to solve 
the Dirac equation with the scalar potential inferred from 
soliton theory. Furthermore, summation of the scalar 
condensates of all continuum states in the Dirac sea and 
the partially filled bound states must be performed to 
check self-consistency. Since the n soliton solutions are 
rather complicated, this might seem hopeless. Remark- 
ably, as we shall show in this section, soliton theory pro- 
vides us with exactly the information needed to perform 
this task in closed analytical form. 

The TDHF Dirac spinor tp for any type I solution sat- 
isfies the extended Dirac equation (|33l34j) . On the other 
hand, the auxiliary spinor tp in the inverse scattering 
problem of the sinh-Gordon equation solves Eqs. (|41I42I) . 
As originally exploited in [10| for the classical N = 1 
GN model and applied to type I solutions of the large N 
GN model in |{|, this implies that the two problems are 
related by a non-Abelian gauge transformation. The lan- 
guage of gauge transformations is adequate here because 
the integrability conditions have the mathematical form 
of a vanishing non-Abelian field strength tensor. Similar 
ideas have been used recently to map the sinh-Gordon 
theory onto string theory in anti de Sitter space AdS3 
fl4| . or the GN model onto string theory Q. We in- 
troduce a gauge transformation CI relating tp and tp as 
follows, 

tp = Clip, 

C x = Cl(U-d)Cl-\ 

C 2 = CliV-d)^ 1 . (49) 

Upon identifying u with In S 2 , we find (modulo an arbi- 
trary normalization factor) 



CI = 



2( 2C 
S -S 



(50) 



With the TDHF spinors at hand, we are now in a position 
to address the issue of self-consistency. Let us start with 



the continuum spinors. Using the gauge transformation 
(|49l50p . wc first write 

^ 1= M2Cp+, il>2=MStp-, (51) 

with tp± = tpx ± ip%. Notice that the linear combinations 



tp± = ±1 - A*(C) 



1 



1 ± a 



A | e i(c *-*> 



(52) 



are actually simpler than 951,2 ■ The normalization factor 
Af will be chosen such as to recover the free Dirac spinor 
d?]) at x — > —00. Using 



lim ipi=Q, lim p 2 = e l ^ z &\ 



this yields 



N = 



(53) 



(54) 



The scalar density can now easily be evaluated with the 
help of Eq. ((151). 



ipip 



4C 



1 + 4c, 2 



:S. 



(55) 



Owing to the vacuum gap equation, the self-consistency 
condition is fulfilled by the negative energy continuum 
states alone, see Eqs. (|11I12[) . It remains to be shown 
that the bound states do not destroy this result. If the 
solitons are far apart, each of them possesses a normaliz- 
able bound state. One therefore expects the presence of 
n bound states in the n baryon problem. Following an 
observation made in Sec. [IT] in the one soliton case, we try 
to generate the bound state spinors from the continuum 
spinors by analytical continuation to imaginary spectral 
parameters. We find that the bound state originating 
from the fc-th soliton can indeed be obtained by setting 

c = a, 



A t/-\ a kl 

MO -> -j-, 
1 



<p- 



1 



-x 



vCT) \ 1 + a )h 
1 < 1 A 



(56) 



The fact that the ±l-terms in (|52|) have disappeared is 
instrumental for the normalizability of the bound states. 
For x — > — 00, A vanishes so that the spinor also vanishes. 
For x — > +00, A increases exponentially but a behaves 
as A 2 , so that again the spinor vanishes. According to 
Eqs. (|39| and (|47| . (k and Cfe(0) have the phase (— i), Xk 
has the phase \f--~i and a k i is real. This shows already 
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that ip + and if- arc in phase. 
Dirac-HF spinor for £ = ( k , 



The components of the 



(57) 



then differ by a phase i so that the scalar density indeed 
vanishes for the bound states. Hence the situation is the 
same as for the single kink: The valence fcrmions play no 
role for the issue of self-consistency The explicit spinors 
■0( fc ) will be needed nevertheless to evaluate the fermion 
density. The only missing piece is the normalization con- 
stant A/fe, to be determined from the integral over the 
density 



dx 



(l4 



0)|2 



1^1 



= 1. 



(58) 



It can easily be found by considering times when the 
solitons are well separated, where it reduces to the one- 
soliton case, cf. Eq. ([15]). 



A4 



(59) 



This completes the proof that the n soliton solution 
of the sinh-Gordon equation yields a self-consistent solu- 
tion of the TDHF equation in the GN model. It covers 
the kink baryon reviewed in Sec. HH and the kink-antikink 
scattering solution of @ as special cases. The n baryon 
solution describes the general scattering problem of an al- 
ternating succession of n kinks and antikinks. Each one 
can carry an arbitrary number of fermions in the allowed 
range and has arbitrary initial positions and velocities, 
parametrized by the constants a k ,v k . The fact that this 
problem can still be solved in closed analytical form, in- 
cluding the polarization of the Dirac sea, is remarkable. 
In the remaining sections we will first cast the results in 
a form more suitable for applications and then discuss 
the physics of the n baryon collision in more detail. 



V. USEFUL EXPRESSIONS FOR SCALAR 
POTENTIAL, SPINORS AND DENSITY 

The preceding section contains all the ingredients 
needed for the full TDHF solution of n interacting kinks 
and antikinks. Yet these results are not yet in a form 
well suited for practical computations with computer al- 
gebra. If one tries to evaluate them, for example with 
Maple, one notices that the number of terms increases 
rapidly with n and algebraic manipulations become pro- 
hibitive for rather small n values already. The aim of 
the present section is to present an alternative formula- 
tion which has proven more convenient for applications. 
It is adapted from a work of Bowtell and Stuart on the 
sine-Gordon equation [3C} and makes the structure of the 
n-soliton solution more transparent. It also facilitates the 
computations of time delays in Sec. I VII and has proven 
to be a prerequisite for practical calculations of sizeable 



number of solitons to be discussed in Sec. IVIII Besides 
developing this approach for both scalar potential and 
TDHF spinors in general case, we have also included in 
this section the proof that the total fermion density is 
proportional to the bound state contribution, generaliz- 
ing Eq. (f26|) to n baryons. This will also be of great help 
for the computations described in Sec. IVIII 

We start with the construction of the n-soliton poten- 
tial S. Since S and — S differ only by a 75 transformation, 
they describe the same physics and we can choose 



lim S = 1 



(60) 



without loss of generality. The single antikink can been 
written in the form 



S = 



l-n 

l + n' 



(61) 



see Eq. (| 14[) . Following the ap pro ach of Bowtell and Stu- 
art in the sine- Gordon case [30J, we first note that n 
non-interacting solitons are described by simply taking 
the product of n onc-soliton solutions, 



fe=i 



l-Tfe 



with 



T fe = exp {27 fe (x - v k t) + 2a k } 



(62) 



(63) 



Clearly, this 2n-parameter ansatz will solve the sinh- 
Gordon equation as long as all solitons are far apart. 
Physically it may be thought of as initial or final con- 
figuration of an n baryon scattering process. S exhibits 
an alternating sequence of n kinks and antikinks. Its 
behavior at spatial asymptotics for fixed time is 



lim S=l, lim S = (-l) r 



X — > — X 



(64) 



Next, we expand the numerator and denominator of 
S. To explain the general construction of the interacting 
soliton solution, it is sufficient to consider n = 3, 



S = 



1 - n - r 2 - t 3 + TiT 2 + T1T3 + t 2 t 3 - nT 2 T 3 

1 + T\ + T 2 + T 3 + TlT-2 + T\Tj, + T 2 T 3 + T1T2T3 ' 



(65) 

Numerator and denominator are multivariate polynomi- 
als of order n in the r k . In order to arrive at the interact- 
ing soliton solution, inspect each monomial of numerator 
and denominator. If it contains T k and tj, multiply it by 
where v k i is the relative velocity of solitons k and I 
(more precisely the velocity of soliton k in the center-of- 
velocity frame of solitons k and I) 



v k i 



1 - v h vi - y/(l-v£)(l-vf) 



Vk - vi 

In our example (n = 3), this prescription yields 



(66) 
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S = 



1 - Tl - T 2 - T 3 + V\ 2 T\T 2 



v\ 3 T 2 T 3 



(v 12 V 13 V 2 3) 2 TlT 2 T 3 



1 + Tl + T 2 + T 3 + v{ 2 TlT 2 + v{ 3 TlT 3 + V^T 2 T 3 + (vi 2 V 13 V 23 ) 2 TlT 2 T 3 



(67) 



r 



This is already the full 3-soliton solution. In other words, 
lnS* 2 solves the sinh- Gordon equation for all values of 
(x,t). Notice that the velocities Vi all have to be chosen 



differently. If Vi 



the result collapses to the n — 1 



soliton solution. More generally, we write the n soliton 
scalar potential as 



(68) 



with 



a¥(t) = l± n , 

A { 1\t) = i±(n + T 2 ) + v 2 12 t iT2 , (69) 

1 ± (Tl + T 2 + T 3 ) + V\ 2 T\T<1 + v\ z TlT 3 
+V 23 T 2 T 3 ± (wi 2 fl3W23) 2 nT2T 3 , 



4 3) (r) 



etc. The relationship between this notation and the one 
in previous sections is made by the following useful equa- 
tions 



Tk 



Ofci 



Ci -c* 



3 



+ Cfc 

= exp 



|2i (a* 



z 



2a fc 



(k + (l 



(70) 



One can now check that the functions A± can equiva- 
lently be expressed as determinants, 



A ( ±\t) = det (l±a< n >) 



(71) 



where aW is the matrix a of Eq. ([70]) for the n soliton 
case. Thus we recover the result (|4"6")l , confirming that 
ln[S , (")] 2 with from Eq. (JBSJI is the n soliton solu- 
tion of the sinh-Gordon equation. The advantage of the 
present algorithm is the fact that it is very easy to im- 
plement in computer algebra and makes the structure of 
the potential more transparent. 

A similar procedure works for the TDHF spinors as 
well. The continuum spinors for the n soliton problem 
can be written as 



1 



2(A { ^( K ) 

(»)/ 



with 



v/TWl-^+V) J a { ; i \t) 



Ci-C 



Ci + C 



(72) 



(73) 



To get the bound state which belongs to the k-th soliton, 
replace the normalization factor in (|72|) by {59} and £ by 



(fc) 



2(A^\k) 



V2iCk V -A^'(k) / ^W( r ) 



(74) 



C=Cfe 



This is significantly simpler than the continuum state, 
since all monomials in the numerators containing a factor 
Kk vanish. 

If one evaluates the fermion densities with Maple for 
small n values, one finds that the simple relation between 
induced and valence fermion density found in the one- 
and two-soliton cases generalizes to n solitons (p^ — 



<p> = £ 



fc=i 



(75) 



Hence one can reconstruct the full fermion density from 
the discrete states alone. Eq. (|75|) can be proven for 
general n with the help of Cauchy's theorem. Since the 
analytic structure of the fermion density is rather compli- 
cated, we demonstrate the corresponding relation for the 
simpler case of the pscudoscalar condensate. As pointed 
out in Sec. HJ the divergence of the vector and axial 
currents establishes a close relationship between various 
fermion bilinears. Eliminating the axial density (ps) from 
Eqs. (|31[) , we can express the fermion density directly in 
terms of the pscudoscalar condensate, 



(76) 



d^(p) = -^25(^75^), 



so that it is sufficient to prove the analogue of Eq. ([75)1 
for (ipij^ip). The pscudoscalar density for a single orbit 
reads 



For a continuum state [see Eq. (|72|)]. we get 

2< 



tpij 5 tp 



1 + 4C 2 



c.c. 



Note the useful relations 
(A^iK)^* =A ! l\k*), 



Q-C 



(77) 



(78) 



(79) 



Wc perform the sum over modes as an integral over the 
spectral parameter, using the residue theorem in the com- 
plex £ plane. The integrand is an analytic, even function 
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of £ falling off like 1 /C 2 at infinity, so that we can extend 
the g?C integration from — oo to +00 and apply Cauchy's 
theorem. In the lower half-plane there are simple poles 
at £ = Cfc arising from ft*,. Since A_ (ft*) is linear in 
each ft*,, we can evaluate the fc-th residue by setting 



A^\k*) 



* d 



XU = 



( "V) 

(«-!) f\* 

l A l,fc) •■ 



= ft£.4 

(the argument A*. fc is missing in A_ 



Ki,k) 



(80) 



When applying 



the residue theorem, A 



(n-l) 



in this expression has to be 



evaluated at the pole £ = Cfe, 



(c, 














-cJ 



so that 



— T„ (81) 



(82) 



This amounts to substituting 



,4^ (ft*) 



7* 



.4L n) (ft) 



(83) 



Inserting the non-singular factors evaluated at the pole 
and summing over all poles at C = Cfcs the residue theo- 
rem gives the following contribution from the continuum 
states to the pscudoscalar condensate 



A^\k)A^\k) 



Tk- 



fe=i 



4V) 



(84) 



For the bound states on the other hand, a straightforward 
evaluation of the pseudoscalar condensate using the wave 
functions (f74|) yields 



bound 



2$> 

fe=i 



^ n) (ft)^i n) (ft) 



fcTfc- 



c=c* 



^ l) (r) 



so that the total condensate becomes 



(#75^) = ( N * ~ Y ) (^7s^) (fe) 
fc=l ^ ' 



(85) 



(86) 



where N^ipi^ip)^ denotes the fc-th term in the sum of 
Eq. ([55]). Due to Eq. (JTBJ this also proves (|75|) . 

Summarizing this section, we note that the most im- 
portant results arc Eqs. (|68|) for the scalar mean field, 
(|72p for the continuum spinors and (|74p for the bound 
state spinors, together with the constructive algorithm 
illustrated in Eqs. (|69[) and the final expression for the 
fcrmion density, Eq. ([75]) . This is the basis for all the 
concrete applications discussed in Sec. IVIII Moreover, 
Eq. (p38"|) is helpful for deriving the asymptotics in Sec. lVIl 



VI. ASYMPTOTICS FOR t -> ±00, PHASE 
SHIFTS AND TIME DELAYS 

The only observable in an clastic scattering process in 
1+1 dimensions is the time delay. Here we shall compute 
the time delay experienced by each of the n baryons. In 
order to determine the time delay, we need the asymp- 
totics of the scalar potential S for t — > ±00. As a 
byproduct, this will teach us how to translate the pa- 
rameters (ak,Vk) into initial positions and velocities of 
the baryons. 

We order the solitons according to the velocities Uj, 

vi > v 2 > ... > v n . (87) 

Then for t — > —00, S describes n incoming (anti-)solitons 
with the functional form (only valid in the vicinity of the 
corresponding soliton) 

l-n 



si 1 ' = 



r(fe) 



i + V 

1 1 \fc+i \ ~ r fc rii=i y ik 

1 + n Ui=i v tk 



(88) 



(fc = 2, n). They are ordered from left to right, start- 
ing with an anti-soliton. For t — > 00, S describes n out- 
going (anti-)solitons with the functional form (again only 
valid in the vicinity of each soliton) 

(89) 



,(i) 



S. 



(fc) 



1 

(-1) 



1 TT n 2 

fc+l 1 ~ T n+l-k llj= n +2-fc v n+l-k,i 

1 + Tn+l-k lli =n +2-k V n+l-k,i 



(fc = 2, ...,n). They are also ordered from left to right, 
starting with an anti-soliton. If one follows the baryon 
density, one finds that it is exchanged in each two-body 
collision. This is a direct consequence of the fact that 
the scalar potential is transparent. Hence a particular 
fcrmion cluster gets transferred from the incoming soliton 
fc to the outgoing soliton n , + 1 — fc; the spatial order 
is inverted. Physically relevant is presumably only the 
time delay for the fcrmion clusters, not the (anti-)kinks. 
This is equivalent to computing the time delay from the 
asymptotic form of S, comparing kinks with the same 
at t — > ±00. The result for fc = 2, n — 1 is 



(At) fc = 



i" (n;u 



In 



(nti 1 ^ 



For fc 



27fc^fc 

1 and k = n, one finds 



(90) 



(At)i = m(nr = x) {At)n 

271 Vl 



(91) 

In the special case of two solitons in the center-of- velocity 
frame, we recover the result of 0] 



v, v 2 



-V, V12 



(At)! = (At)a 



lnu 2 

27U 



(92) 
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In the soliton literature, one also introduces a "phase 
shift" related to the time delay by [3l[ 



6 k = -2"f k v k (At) k . 



(93) 



The total phase shift is the sum of the phase shifts in- 
duced by independent collisions with all other solitons. 
The 6 k satisfy 



(94) 



k=l 



To specify the initial conditions, it is helpful to note 
the equation of motion of the fc-th incoming soliton, 



a k + In C k 
27fc 



x = v k t - 
C l = 1, 

k-l 

Ck = ]Jvl, k = 2,...,n. 



(95) 



Similarly, the equation of motion of the fc-th outgoing 
soliton (numbered in inverse order, i.e., according to the 
fcrmion clusters they carry) reads 



V k t 



a k +\nC' k 
27fc 



C„ — 1 , 



Cfc = J| Wfci, k = 1, ...,n- 1. 



(96) 



Denoting the initial time by t = —T, the initial positions 
of the solitons are given by 



(fc) m a/t+lnCfc 

x =-VkT — . 

^7fe 



(97) 



This tells us how to choose the parameters , given the 
initial soliton velocities and positions, namely as 



a k = -lnC fc -2 7fc (4 fe) + v k T^j 
with C k from Eq. ([95]). 



(98) 



VII. ANATOMY OF THE N BARYON 
COLLISION AND ILLUSTRATIVE EXAMPLES 

What happens if one prepares n alternating, well sepa- 
rated kinks and antikinks with different initial velocities, 
carrying different numbers of fermions or antifcrmions? 
We are now in a position to predict the time evolution of 
this initial configuration in the GN model. In general, it 
would be very hard to characterize such a complex col- 
lision process. In our case there are several simplifying 
features which enable us to draw a full picture. 

The initial and final states of an n-body collision may 
be described as in the previous section — the solitons are 




FIG. 1: Schematic drawing of generic multi-soliton collision 
for n — 4. Time t runs vertically, the x-axis is horizontal. 
Fermions travel approximately along the straight lines, in- 
tersection points denote two-soliton collisions. Every soliton 
scatters exactly once from every other soliton. 




-10 



-20 



FIG. 2: Time evolution of scalar mean field S for the 
4 soliton case sketched in Fig. [1] Parameters: a — 
{50.5, 25.2, 21.8, 14.1}, v = {0.9, 0.7, 0.4, 0.1}. 



P o 




FIG. 3: Like Fig. [51 but fermion density shown. Solitons 1 
and 4 have maximal fermion number N/2, solitons 2 and 3 
are empty. 
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widely spaced and ordered according to their velocities, 
the fastest one being leftmost in the incoming and right- 
most in the outgoing state. We illustrate such a process 
schematically in Fig. [1] for the case of n — 4 solitons, 
in a frame where all velocities Vi are positive. Since a 
kink and an antikink cannot pass through each other, it 
looks as if the solitons repel and stay in the same order. 
However, due to the fact that the self-consistent poten- 
tial is transparent, the fcrmions carried by each kink or 
antikink can only move forward. In every two-soliton 
collision, the fermions get exchanged as discussed in 
Inelastic processes are suppressed due to the integrability 
of the GN model. In Fig. [T] the fermions move roughly 
along the straight lines (ignoring interaction effects). The 
intersection points of two straight lines signal two-body 
collisions. Obviously, every baryon interacts with every 
other one exactly once. The complete time evolution of S 
including interaction effects is shown in Fig. [3J where one 
recognizes time delays. Fig. [3] shows the corresponding 
time evolution of the fermion density. To simplify the pic- 
ture, we have assumed that solitons 1 and 4 have maximal 
fermion number N/2, whereas solitons 2 and 3 carry no 
fermions at all. We see that the fast, Lorentz contracted 
fermion cluster of soliton 1 passes through the collision 
zone almost unaffected. The wider peak corresponding to 
the slower fermions of soliton 4 suffers stronger interac- 
tion effects, being also scattered by the "empty" solitons 
2 and 3. If we had loaded any of the solitons with an- 
tifermions by choosing an occupation of the valence level 
< N/2, we would observe that fcrmions and antifcrmions 
also pass through each other, due to the absence of an- 
nihilation processes. Note also that the graph shown in 
Fig. [2] is independent of the fermion content of the soli- 
tons. It would even hold in the case where all solitons 
have vanishing fermion number, so that neither baryons 
nor bosons are involved. Nevertheless we would be deal- 
ing with a valid solution of a quantum field theory. This 
underlines the non-perturbative character of the whole 
approach. 

Let us now consider some further illustrative exam- 
ples. Our original motivation for studying the GN 
model came form strong interaction physics. In real 
life, natural many-baryon problems would involve nu- 
clei. If the 2-soliton scattering is taken as a toy model 
for nucleon-nucleon scattering, one would like to ad- 
dress next nucleon-nucleus or nucleus-nucleus collisions 
at the elementary fermion level. Unfortunately, the GN 
model has no "nuclei", i.e., bound states of baryons. The 
baryon-baryon interaction is repulsive. "Nuclear matter" 
exists in the form of a soliton crystal, but it is neither self- 
bound, nor does it saturate. Scattering problems with 
more than two incident particles on the other hand have 
no obvious analogue in particle physics. Therefore, the 
best we can do to mock up nuclear targets or projectiles 
in our toy world is to use "trains" of solitons with nearly 
equal velocities. Although unstable, such a configuration 
will stay together for a time long enough to study scatter- 
ing processes. These trains of solitons may be thought 




FIG. 4: Relativistic "proton-nucleus" collision simulated by 
the collision of a single soliton with a train of 5 solitons, ap- 
proximatively at rest ("labarotory frame"). Time evolution of 
scalar potential is shown. Solitons behave like hard spheres. 
Parameters: a, = 0, v = {0.9, 0.02, 0.01, 0, -0.01, -0.02}. 




FIG. 5: Like Fig. [4] but fermion density shown. All projec- 
tile and target solitons carry the allowed maximum of N/2 
fermions. 



of as chunks of soliton crystals ("nuclear matter"). In 
applications of the present model to other fields like con- 
densed matter physics, the interest would presumably be 
in a different kind of n-soliton problem. The formulae 
given in Sec.|V]should enable the reader to produce easily 
any desired result by choosing appropriate parameters. 

Proceeding in this spirit, we show in Figs. 2] and [5] 
an example of the analogue of a baryon-nucleus collision 
for 1+5 solitons, in the (approximate) rest frame of the 
target "nucleus". The kinks behave much like classical 
hard spheres, i.e., the incoming projectile gets stopped 
when it hits the first target baryon, and the last target 
baryon leaves, carrying away the momentum. This can 
be inferred from the scalar potential in Fig. 2J To illus- 
trate the fate of the fermions, we fully load the projectile 
and target baryons with N/2 fermions each. As shown 
in Fig. [5j the fast projectile fermions then hop from one 
soliton to the next one repeatedly during the collision, 
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FIG. 6: Relativistic "nucleus-nucleus" collision simulated by 
the collision of two trains of 4 solitons each, in the center- 
of-mass frame. Scalar potential shown. Parameters: a.% = 0, 
v = {0.95, 0.94, 0.93, 0.92, -0.92, -0.93, -0.94, -0.95}. 




FIG. 7: Like Fig. [6] but fermion density shown. Each soliton 
carries N/2 fermions. 



until they emerge in the emitted, rightmost soliton and 
move along with it. 

Owing to the relativistic invariance of the formalism we 
can study these collision processes in any desired Lorentz 
frame. In our last example, we choose the center-of-mass 
frame of two "nuclei", each one consisting of 4 solitons 
carrying the maximal number of fermions. This is the 
closest we can come to simulate a "relativistic nucleus- 
nucleus collision" in the GN model. Figs. [5] and [7] show 
again that the solitons repel each other, whereas the 
fermions keep moving forward. A combination of inte- 
grability and transparent mean field is behind this simple 
scenario. 

By comparing Fig. |6] for the scalar potential and Fig. [7] 
for the density, one may be tempted to conclude that 
both figures show multi-soliton collisions. Indeed, in both 
cases all the lumps emerge unchanged from the collision 
process. However this interpretation is only valid for S 
and the underlying sinh-Gordon equation. That the den- 
sity has no solitonic character already follows from the 



fact that the normalization of each fermion cluster can 
be arbitrarily chosen. Formally, whereas S is the solu- 
tion of a non-linear differential equation, the density can 
be thought to arise from a linear equation where S en- 
ters as an external field, similar to the spinors in the 
TDHF equation. In any case, the fact that our solitons 
carry fermions is an interesting aspect not shared by stan- 
dard applications of solitons in physics, but reminiscent 
of early soliton bag models [32j in 3+1 dimensions. 



VIII. SUMMARY AND OUTLOOK 

The observation that solitons share some properties 
with elementary particles is as old as soliton theory. In 
the GN model, this relationship can now be made very 
precise. The underlying quantum field theory is purely 
fermionic. It produces dynamically multi-fcrmion bound 
states. In the large N limit, the appropriate scmiclassi- 
cal setting is the relativistic HF approach. The scalar HF 
potential is a classical field with solitonic character, but 
the bound fermions are also relevant for understanding 
the structure of "hadrons" . This is of course well known 
since long time already. The new insight which we can 
add now is the fact that for a certain class of particularly 
simple solutions (called type I), the whole dynamics can 
be decoupled from the fermions and cast into the form 
of a non-linear differential equation for the scalar mean 
field. This equation turns out to be the sinh-Gordon 
equation. Apparently one can bypass the complicated 
self-consistency issue for these particular solutions and 
arrive at the self-consistent solution by just solving a sin- 
gle, nonlinear differential equation for the "master field" 
S. The fermions then follow the motion of the solitons, 
but do not react back in any way. Since the relevant 
soliton equation is well known, this enabled us to solve a 
rather involved problem in closed analytical form, namely 
the dynamics of n kink and antikink baryons with arbi- 
trary fermion number, initial positions and velocities. We 
have analyzed this type of scattering process and found 
that it has many unrealistic features from the point of 
view of particle physics. However, here we have no choice 
since we are not dealing with a phenomenological model, 
but the solution of a given quantum field theory, Eq. (QJ , 
in the large N limit. Actually, examples in theoretical 
physics where the dynamics of a number of composite 
particles can be analyzed exactly at the elementary con- 
stituent level are extremely rare, even in non-rclativistic 
many-body physics. In our case, Lorentz covariance is 
strictly maintained and the polarization of the Dirac sea 
fully taken into account. 

The methods developed here in a concrete example 
may have some potential for generalizations. One strik- 
ing observation is the fact that the TDHF spinors are 
apparently closely related to auxiliary spinors introduced 
in soliton theory when one looks for solutions via the in- 
verse scattering method. It is very likely that there is a 
more general principle behind this apparent coincidence. 
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It was certainly important that we restricted ourselves 
to type I solutions of the TDHF equations. All other 
analytically known solutions of the massless or massive 
GN model are actually type II and therefore related to 
the N = 2 classical GN model. It would be interest- 
ing to generalize our approach to this more general case, 
thereby extending the pool of exact solutions, perhaps 
even to non-intcgrable field theories like the massive GN 



model. 
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